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Abstract. We have investigated some of the properties of dense sub-nuclear matter at the crustal region 
(both the outer crust and the inner crust region) of a magnetar. The relativistic version of Thomas- 
Fermi (TF) model is used in presence of strong quantizing magnetic field for the outer crust matter. The 
compressed matter in the outer crust, which is a crystal of metallic iron, is replaced by a regular array of 
spherically symmetric Wigner-Seitz (WS) cells. In the inner crust region, a mixture of iron and heavier 
neutron rich nuclei along with electrons and free neutrons has been considered. Conventional Harrison- 
Wheeler (HW) and Bethe-Baym-Pethick (BBP) equation of states are used for the nuclear mass formula. 
A lot of significant changes in the characteristic properties of dense crustal matter, both at the outer crust 
and the inner crust, have been observed. 
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1 Introduction 

Magnetars are the most exotic stellar objects, believed to 
be strongly magnetized young neutron stars. The surface 
magnetic field for such objects are observed to be > lO^^G 
[UllllaJll]. Then it is quite likely that the field at the inte- 
rior, even at the inner crust region may be stronger than 
the surface value (can be predicted theoretically by scalar 
Virial theorem) . If the internal field strength is happened 
to be so high, then most of the physical and chemical 
properties of dense stellar matter of the magnetars must 
change significantly from the conventional picture [SUHIIZ] • 
A lot of investigations have already been done on the ef- 
fect of strong quantizing magnetic field on various physical 
properties of dense stellar matter inside neutron stars as 
well as quark matter inside quark stars or hybrid stars, 
including the efi^ect on quark-hadron phase transition at 
the core region of a compact neutron star [S] . The effect of 
such strong magnetic field on various elementary processes 
inside neutron stars and quark stars or hybrid stars have 
also been studied. The effect of strong quantizing mag- 
netic field on the /3-equilibration among the constituents 
have also been investigated [9,. It has also been shown 
that strong quantizing magnetic field acts like a catalyst to 
generate fermion mass dynamically, i.e., chiral symmetry 
breaking occurs in presence of strong quantizing magnetic 

field uoiimiii]. 

In this article we present our investigation on the effect 
of strong magnetic field on the crustal matter of magne- 
tars. The work is divided into two parts: in the first part. 



based on one of our very recent work |13j , we have inves- 
tigated the effect of strong quantizing magnetic field on 
the outer crust matter. In the second part, we have stud- 
ied the properties of compact sub-nuclear matter at the 
inner crust region in presence of such strong quantizing 
magnetic field. 



The paper is organized in the following manner: In 
section 2, the effect of strong quantizing magnetic field on 
the crustal matter of magnetars is discussed. In the outer 
crust region, the matter with dense crystalline structure 
of metallic iron at sub-nuclear density are replaced by an 
array of spherically symmetric WS cells with positively 
charged nuclei at the centre surrounded by non-uniform 
dense electron gas with over all charge neutrality. In the 
inner crust we have used the conventional HW and BBP 
equation of states for the nuclear matter [HIITSIITC] , con- 
sisting of iron along with some more heavier neutron rich 
nuclei. In one of our forth coming article we shall report 
the effect strong magnetic field on inner crust matter us- 
ing some modern type nuclear mass formulas and make a 
comparative study with those models. Here we have pre- 
sented with detailed numerical computation the effect of 
strong magnetic field on the inner crust matter, which as 
mentioned above is assumed to be a mixture of iron, some 
heavier neutron rich nuclei, which we found to be spe- 
cially true in presence of strong quantizing magnetic field, 
electrons and free neutrons. The presence of free neutrons 
are considered beyond neutron drip density. Finally, in the 
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last section, we have given the conclusions and discussed 
the importance and future prospects of the present work. 



2 Crustal Matter 

In the first part of this section, we present our study on 
the effect of strong magnetic field on the outer crust region 
with the composition mentioned in the introduction part. 
For a typical neutron star, the width of the outer crust 
region is ^ 0.2 — 0.4Km. In a recent work we have devel- 
oped an exact formalism for the relativistic version of TF 
model in presence of strong quantizing magnetic field [13| . 
This formalism is used within the limitation of TF model 
[T7] , to obtain the equation of state for crustal matter of a 
typical magnetar. In this model it has been assumed that 
for a WS cell, the atomic number is Z and the mass num- 
ber is A. To make each cell electrically charge neutral, Z 
must also be the number of electrons inside the WS cells. 
In this model, the modified form of Poisson's equation 
is given by [13, 






1^=0 



{2-5,o){<t>\x)~^lx^Y'^ 



(1) 



where (t){x) is the modified form of electrostatic poten- 
tial related to the original Coulomb potential V{x) by the 
relation 

eV{x)l (2) 



4){x) = Y^W 



with ^e, the electron chemical potential, assumed to be 
constant throughout the cell (this is the so called Thomas- 
Fermi condition), the dimensionless scaled radial coordi- 
nate X is related to the actual radial coordinate r of WS 
cells by the relation: x = r//i, with 



A* = 



2e^B 



1/2 



(3) 



e — [e|, the magnitude of electron charge, B is the constant 
magnetic field, assumed to be along z-direction (we have 
chosen the gauge A'' = (0, 0, xB, 0)), 



~Z^' 



(4) 



with rrii, — (m^ + 2zye-B)^/^, mg — O.SMeV, the electron 
rest mass and z^ = 0, 1, 2, , i^max, is the Landau quan- 
tum number for the electrons, I'max is the upper limit of 
Landau quantum number (the upper limit of the Landau 
quantum number i^max is finite at zero temperature, other- 
wise fc'max — oo). Finally, the factor {2 — 5uq) indicates that 
the zeroth Landau level is singly degenerate, whereas, all 
other states are doubly degenerate. In this article we have 
assumed that the electron gas in the dense crustal mat- 
ter of crystalline metallic iron is strongly degenerate and 
is considered to be at zero temperature. Then it is quite 
obvious from the non- negative value of electron Fermi mo- 
mentum, that the upper limit of Landau quantum number 
t-max is given by 



^^e 



2eB 



(5) 



To obtain numerical solution for (pix) for a given mag- 
netic field and a particular set of Z and A, we further 
assume that instead of a point object, the nucleus at the 
centre of a WS cell, has a finite dimension and is assumed 
to be spherical in nature, so that the corresponding radius 
r„ = tqA^^^, with To — 1.12fm. Such a choice also removes 
the singularity problem of TF equation at the origin [19| . 
Therefore it is not necessary to implement the prescrip- 
tion given by Feynman, Metropolis and Teller to obtain 
the numerical solution for TF equation [20 . Again from 
the physics point of view, the potential must satisfy the 
boundary conditions, given by 



rV{r) = Ze for 



and 



dV_ 
dr 



for 



(6) 

where Vg is the radius of the WS cell. Then by simple 
algebraic manipulation it is easy to show that 4'{x), the 
modified form of Coulomb potential satisfies the boundary 
conditions 



0(a;)|x=x„ 



1 and -— I 
dx 



cj>{x) 



(7) 



where Xn = ?'n/M, the scaled nuclear radius and Xs = ?'s/a'j 
the corresponding scaled radius of the WS cell. Both these 
quantities are dimensionless. 

Further, the right hand side of the Poisson's equation 
(eqn.(l)) must be real. Which requires 0^a; < |0(cc)|. This 
is an additional condition, to be satisfied by the upper 
limit of Landau quantum number, and may be written by 
the following inequality: 



.{x)< 



-4>{xf 



2eB 



(8) 



From the definition of Landau quantum number, the above 
inequality must necessarily be > 0. The above equation 
also shows that the upper limit of Landau quantum num- 
ber depends on the position (a; or r coordinates) of the 
electron within the WS cell, with which it is associated. 

Since electron distribution is non-uniform within each 
cell, the Fermi momentum of a particular electron must 
depend on its positional coordinate in the cell. Then it is 
expected that the variation of pi?(r) will be such that the 
electron chemical potential, given by 



/^. 



-t- eV{r) = {pF{rf + ml + 2v{r)eBY''^ 



remain constant throughout the cell with proper space 
dependent Landau quantum number v. 

Satisfying all these conditions, of which some of them 
are particularly necessary for this model, we have solved 
the Poisson's equation numerically within the range of r 
from nuclear surface to the WS cell boundary, for B = 

B'f\ 10 X b[^\ 10^ X b[^^ and 10^ x b[^\ where B^^'' is 
the quantum critical limit for the magnetic field at and 
above which the Landau levels are populated for the rel- 
ativistic electrons, and may be expressed in terms of the 
universal constants, given by eBc — ml (with the choice 
of unit h — c — 1). The magnetic field beyond this limit 
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is called the quantizing magnetic field and the quantum 
mechanical effect plays an important role in this domain. 
The numerical value of this field is 4.43 x lO^^G. Within 
the range Xn < x < Xs, the numerical solution for ^(a;) 
can be fitted exactly by a straight line for a particu- 
lar magnetic field strength B and can be expressed as 
4>{x) = ax + b. Obviously, the parameters a and b are 
functions of magnetic field strength. We have obtained 
this functional form by minimizing x^ numerically using 
a computer code for multi-parameter fitting program. In 
Table-I we have shown the variation of the parameters a 
and b with magnetic field strength. 

Because of some numerical in-accuracy, the values of 
the parameter b, the intersections of the straight lines with 
y-axis, are not exactly one. However, for all the magnetic 
field values, the parameter b is very close to one. The nu- 
merical values for the other parameter a is always negative 
and its magnitude increases with the increase in magnetic 
field B. Which actually means that the radius or the vol- 
ume of a particular WS cell decreases with the increase in 
magnetic field strength. However, in this particular article, 
we have avoided the origin, considering a finite size nucleus 
of radius r„ at the centre. Then we have on the nuclear 
surface 6=1 — arn/ ^. In Table-I we have also shown the 
explicit variation oixs, the scaled surface radius of the WS 
cells, for various magnetic field strength. While solving the 
Poisson's equation numerically for a given magnetic field 
strength B, the instruction is given in our numerical code 
to terminate if the surface condition, given by eqn.(7) is 
satisfied and hence we obtain Xg (also Vg = Xg^j) as a func- 
tion of magnetic field strength. In fig.(l) we have shown 
the variation of Xg and also the corresponding Tg as given 
in eqn.(9) below, with the strength of magnetic field. This 
figure clearly shows that the WS cells become more com- 
pressed if the magnetic field becomes stronger. This is in 
some sense analogous to what is called the magnetostric- 
tion in classical magneto-statics. These two curves can also 
be fitted by the power law functions, given by 



Xg = 29.43 X 



rg = 1.704 X 



B 

w 

B 



-0.41 



MeV"^ and 



B. 



{e) 



-0.91 



A 



(9) 



Knowing the scaled Coulomb potential (/'(x) at various x 
points (in the range Xn to Xg) within the WS cell for a 
given magnetic field strength, we have evaluated numeri- 
cally ^'niax(a^) at every x points within the cell. The vari- 
ations of ^'inax(2:) with X for four different magnetic field 
strengths are shown in fig. (2). Although ^'max(a;) must be a 
set of discrete numbers, for the sake of illustration, we have 
plotted it as a continuous variable. In this figure, curves 
a and b are for B — Be and 10 x Be respectively. For 
these two curves, the variables i'max(a;) is plotted along the 
left side y-axis and x (this is actually iix in Mev~^) along 

X-axis at the bottom. Similarly for i? = 50 x Be and 
100 X Be', the variations are shown by the curves a and 
/3 respectively. In this case, i^max(a;) is plotted along the 



right side y-axis and x is plotted along upper x-axis. From 
this figure, it is possible to make a number of conclusions: 
(i) For larger B values, Xg are smaller, (ii) For smaller 
B, Vmax{x) starts with quite large value near the nuclear 
surface, e.g., = 124 and = 15 as shown in curves a and b 
respectively. Whereas, for curve a it starts with I'max — 2 
and for /3, the starting value is fmax = 1- (hi) The discrete 
nature of fmax is obvious from the high magnetic field 
curves a and /?. (iv) Finally, for all the field values, the 
upper limit t'max(a;) becomes exactly zero at the surface 
of the WS cells. In other wards, we can say that all the 
electrons near the WS cell surface are strongly polarized 
and the spins are anti-parallel to the direction of magnetic 
field. This is, of course, a purely relativistic effect. The pos- 
sibility of fully polarized scenario can easily be obtained 
from the analytical solution of Dirac equation for electrons 
in presence of strong quantizing magnetic field. The eigen 
functions will not be simple spinor solutions, whereas the 
energy eigen value will be i?y = (p^ + 2t^ei3 + 7TT-g)^/^, with 
2v = n + \ + TTig, where n = 0, 1, 2, .. is the Landau prin- 
cipal quantum number and nig = ±1, the eigen values for 
the spin operator az pTlfTB] . Hence, for z/ = t'max = 0, 
the only possible choice is the combination n = and 
nig = — 1. Which actually means that in the zeroth Lan- 
dau level the spins of all the electrons are in the direction 
opposite (this is due to negative charge carried by the 
electrons) to the external magnetic field. We have further 
noticed that beyond the field value 100 x i?c (which is 
slightly > lO^^G), the upper limit i^max(a;) becomes iden- 
tically zero not only at the surface region of WS cells, but 
at all the points inside the cell. 

To obtain the density distribution of electrons within 
the WS cells, let us consider the expression for electron 
number density, given by 



eB 



l^max(a;) 



(^) = ^ E (2-MPF(a;) (10) 



where the electron Fermi momentum pf{x) can be ob- 
tained from the TF condition and is given by 



Pf{x) = < Z e 




(11) 



Since the electron density is larger near the nuclear 
surface inside the WS cells compared to the boundary 
values, the effect of electron density dominates over the 
influence of magnetic field strength at the central region, 
whereas near the WS surface, the magnetic field will play 
significant role. This will naturally make the upper limit 
i^max large enough (except for B > 10 '^^G) near the nu- 
clear surface and quite small or identically zero near WS 
cell boundary. The results are also quite obvious because 
of the electron Fermi momentum, which is minimum at 
the nuclear surface and maximum near the the WS cell 
boundary. Later we shall show that the electron kinetic 
energy will also behave like Fermi momentum. Therefore, 
the quantum mechanical effect for low and moderate val- 
ues magnetic field strength will be more significant near 
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the surface region compared to the core region, where the 
effect is almost classical. This is particularly obvious from 
the curves a and b of fig. (2). For these curves, v-may.ix) is 
naturally started with quite large values. 

Using the numerically fitted form of 4){x) we obtain 
Pf{x) for a given magnetic field strength and the cor- 
responding ne{x). In fig. (3) we have plotted the varia- 
tion of electron density within WS cells for four differ- 
ent magnetic field strengths. In this figure, the curves in- 
dicated by the numbers 0,1,2 and 3 are for B = Be , 
B = lOx B'f\ B = Wx B^"^ and B = 10^ x BJf^ respec- 
tively. The physical reason behind such increase in elec- 
tron density is because of the decrease in WS cell volume 
with the increase in magnetic field strength, whereas, the 
total amount charge within the cell remain constant. It is 
quite obvious from these curves that the electron density 
increases with the increase in magnetic field strength. Fur- 
ther, for all the values of B, electron density is maximum 
at the centre and minimum at the surface. The oscillatory 
nature of ne{x) is because of the crossing of various Fermi 
levels within the cell. This is a kind of de Haas van Alphen 
effect. In our model we have assumed eqn.(lO) as given in 
P^ . One must use the condition [35] 



ee(r)-eT/(r)+^ex(r) = He 



constant 



(12) 



then re-cast the differential equation as given by eqn.(l) 
in this article. Here E^x is the electron-electron exchange 
energy per particle. However, with this choice it will be 
an almost impossible task to solve the problem even nu- 
merically. Hence we have made approximation by using 
eqn.(lO) of [H]. 

Let us now evaluate the variation of (i) electron ki- 
netic energy, (ii) electron-nucleus interaction energy, (iii) 
electron-electron direct interaction energy and (iv) electron- 
electron exchange energy within the WS cells. The kinetic 
energy part for electrons at a particular point (r) within 
the WS cell is given by 

r „ i^max(r) 

Eke{x)= / d\— J2 (2-M 



I 



PFir) 



v=0 

,2 I ™2\l/2 



dp,[{pi + mty^'-me] (13) 



Evaluating the integral over p^ analytically and then sub- 
stituting for pf{x) from eqn.(ll), with the numerically 
fitted functional form of 0(a;), one can obtain the elec- 
tron kinetic energy Eke, as a function of r or a; by the 
numerically evaluating the above integral. We have seen 
that within the WS cells, for constant B, the kinetic en- 
ergy part satisfies a power law of the form Eke = cex^ , 
where the parameters a and /? are functions of magnetic 
field strength. In Table-I we have shown these variations. 
The change is more significant for a than /3. Therefore, 
just like the electron density, the electron kinetic energy 
is also an increasing function of radial coordinate x but 
does not show oscillation. 

Next we consider the three possible types of interac- 
tion potential. Let us first consider the electron-nucleus 



interaction part, given by |13j 

Een{r) = -Ze^ / dV— 
= — 47rZe fjL I xdxne{x) 



(14) 



To obtain this quantity within the WS cell for a constant 
B, we substitute the expression for electron number den- 
sity ne{x) from eqn.(lO), and then numerically integrate 
over X In fig. (4) we have shown the variation of —Een{x) 
(which is a positive quantity) as a function of x within 
the WS cell for three different magnetic field strengths: 
B = B^"^ , 10^ X Bi"^ and 10^ x B^"^ . For very low x values, 
i.e., almost on the surface of the nucleus at the centre, the 
magnitude of the potential energy is an increasing func- 
tion of X, which actually means that the region is strongly 
attractive in nature. Again at the surface region, just at 
the skin of the cell, it is again an increasing function of 
X (with larger gradient), so that the attractive field again 
becomes extremely strong at the surface region to keep the 
outer most electrons confined within the cell. At the mid- 
dle region, there is a kind of saturation and force due to 
electron- nucleus interaction almost vanishes. This is quite 
analogous to the normal metallic scenario, in which elec- 
tric field can not exist. 

Next we consider the electron-electron interaction. Let 
us first evaluate the direct term. It is given by [13] 



2 



d rne{r) 



d\' 



ne{r')- 



1 



(15) 



Assuming r as the principal axis and 9 is the angle be- 
tween r and r', we have d'^r — 4tt r^ dr, 
(fir' = 2Trr' "dr' sin 9 d9 (we have assumed that the vec- 
tors r and r' are on the same plane) and |r — r'| = 
(r^ + r' — 2rr' cos 6*)^/^. The limits for both r and r' are 
from r„ to r^ and the range of 9 is from to tt. Then evalu- 
ating the trivial integral over 9 the direct interaction, part 
reduces to the following simple form: 



Eii'> ^Se^TT^l rdrueir) r'"^ dr'n^ir') 

+ / r^drueir) f r'dr'ne{r')\ (16) 



Following the same procedure adopted for E^n case, we 
have evaluated numerically the above coupled integrals. 
In fig. (5) we have shown the variation of E^^{x) within 
the WS cell for three different magnetic field strengths: 
B = B^e\lQ^ X B^^'^ and 10^ x B'f\ Surprisingly, the 
variations are almost identical with \Een\- However, the 
magnitudes are several orders less than the correspond- 
ing electron-nucleus part. Of course, unlike the E^n, the 
direct term is positive throughout and for all the values 
of magnetic field strength. The same kind of qualitative 
nature of | Een \ and E^J (figs.(4)-(5)) makes the system 
of electron gas within the WS cell energetically stable. Of 
course the kinetic energy part and the exchange part also 



Please give a shorter version with: \authorrunning and Xtitlerunning prior to \maketitle 



play important role. The combination of kinetic energy 
and three interaction potentials of which electron-nucleus 
interaction part and the exchange part are negative, gives 
a constant energy for a particular electron at any point 
within the WS cell. This constant energy is the electron 
chemical potential and is also the minimum possible en- 
ergy at a particular point. 

Next we shall consider the exchange term, which is 
negative in nature. The magnitude of exchange energy in- 
tegral corresponding to the ith. electron in the cell is given 
by[I3] 

where the spinor wave function Tp(r) is given by eqns.(2)- 
(5) in [TB] and '0(r) = -0^^)70! the adjoint of the spinor 
and 7o is the zeroth part of the Dirac gamma matrices 7^ . 
Now it is very easy to show that for t = t' 



i;i{r)'ipi{r') 



2m 



LyLzEi/ 



exp[-i{py{y - y') + Pz{z - z')}] 



{Iu:p, {X)lu-py {x') + /,_i;j,^ (x)/._i;p„ {x')} (18) 

Similarly, we have 

— ZiTfi 

'^j{r')i'3{r) = j^j—^ eM^{p'y{y " 2/') +pU^ ^ z')]] 

(19) 



{Iu':p' {x)Iu';p' {x') + Iu'-l:p' {x)I^' ^i-w {x')} 



where I^-^p^ {x) is same as 7^, given by eqn.(5) in [T3]. When 
these two terms are combined, we have, after replacing the 
sum over j by the integrals 



E^" 



LyLz I dp' 



4m"' 



hPF 



dp'z 



PF 



^I^LILIEJ^^ 



(2 - 5,,o) 



j ... j Lydp'yLzdp'zd\d\'^y^^ 

e^v[-i{{Py - P'y){v - V') + [Pz - p'z){z - z')}] 
[{l>^;Py {x)Iu,py {x') + Iu-i,py {x)h-i,p^ {x')} 
{Iu',p'^{x)I^>,p'^{x) + I^'^i,p,^{x)I^'^i,p'^{x')}] (20) 

It is possible to evaluate the integrals over y' and z', given 
by [^ 



+00 



— 00 J —CO 



dy dz 



\r — r 



exp[-i{(py -py,)iy-y') + {pz - Pz'){z - z')}] 



47r 

— exi>{-K\x-x'\) 



(21) 



where K — [{py — Py'Y + {Pz ^ Pz'YY^"^- Further, the 
integral over y and z is given by 



00 /'+00 



dydz ^ LyL 



y-^z 



Since the final form of exchange integral is multidimen- 
sional in nature with an extremely complicated structure 
of integrand, it is absolutely impossible to have analyti- 
cal solution for e — e exchange potential. We have there- 
fore adopted the multi-dimensional Monte-Carlo integra- 
tion code to evaluate the exchange term. It is found that 
the variation of the magnitude of exchange energy with 
the scaled radius x within the WS cells for a constant B 
can be expressed by the power law formula, given by 

\Eir^\=px'>, 

where the parameters p and q are again functions of mag- 
netic field strength B. The variation of the parameters p 
and q with magnetic field are shown in Table-I. 

From the variation of the parameters p and q with B, it 
is quite clear that the magnitude of exchange energy also 
increases with x and also with the magnetic field strength 
B, i.e., minimum at the central nuclear surface and max- 
imum at the WS cell boundary. 

The overall pressure term can also be obtained for the 
electron gas from the total energy Etot, given by 



Pix)=nl{x) 



dEu 
drie 



Pk{x)+Pc{x) 



(22) 



where 



is the electron kinetic pressure part. The momentum inte- 
gral for the kinetic pressure can very easily be obtained, 
and is given by 



Pk 



eB 

2^ 



u=0 

m' In 



^(2-(5,o) PF{pl + ml) 



2\l/2 



2, lPF + {pl + nilf'^ 



"u 



(24) 



Since pp — Pf{x), this expression will give the electron 
kinetic pressure at various points within the WS cell. The 
other term, Pc, coming from the interaction part, E^ — 



E. 



Eee + E^T' , and is given by 

Pe(x)=r.^(^)|^, 



(25) 



00 J —00 



The interaction part of electron pressure is obtained nu- 
merically by substituting the expression for pf{x) and 
evaluating the derivatives over rie numerically for a given 
i3, at a particular point x within the cell and also for a 
given set of A and Z. We have noticed that at the central 
region of WS cells, since the density dominates over the 
magnetic field part, the total pressure is positive, while 
at the outer end, the surface region, since magnetic field 
dominates over the density contribution, it becomes nega- 
tive. It is found numerically that at some point within the 
WS cell for a given B, the total pressure becomes exactly 
zero. It is found that this critical position is a function of 
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magnetic field strength and decreases with the increase in 
magnetic field strength. In fig. (6) we have plotted this vari- 
ation. The functional dependence can also be expressed as 



Xc = a + P exp 



B 



-7- 



B. 



(e) 



MeV~^ 



(26) 



with a = 0.38, /? = 1.28 and 7 = 0.171. If the quantum 
mechanical effect of strong magnetic field dominates over 
the density effect, electrons either occupy only the zeroth 
Landau level or low lying levels. Which further causes a 
decrease in pressure term. As the magnetic field becomes 
strong enough, the density effect vanishes, even near the 
centre, which makes Xc extremely small. 

In the second part of this section we shall study the 
effect of strong quantizing magnetic field on inner crust 
matter at sub-nuclear density. In the conventional neutron 
star model, for a typical neutron star of radius « 10km, 
the width of inner crust is about 0.5 — 0.6km. It is also 
believed that the free neutrons in the inner crust region 
may be in super-fiuid state. The density range of this zone 
depends on the type of equation of state or the mass for- 
mula for the heavy nuclei, including the highly neutron 
rich nuclei. The density range in the conventional picture 
is ~ 10^ — lO^^g/cm^. To investigate the properties of in- 
ner crust matter, here also we replace the metallic iron 
atoms and other heavier atoms with neutron rich nuclei, 
by WS cells as has been done for the outer crust matter. 
We further assume that the free neutrons in this region, 
above the neutron drip point are in normal fiuid state. 

In this part, to study the effect of strong quantizing 
magnetic field on inner crust matter, we consider two types 
of conventional nuclear mass formulae, which are gener- 
ally used below the nuclear saturation density, but appli- 
cable near neutron drip region. These two mass formulae 
are (i) Harrison- Wheeler (HW) equation of state, and (ii) 
Bethe-Baym-Pethick (BBP) equation of state [nill5lll6| . 
We have already mentioned that in a future communica- 
tion the effect of strong magnetic field on the equation of 
states with more modern type nuclear mass formulas will 
be reported. 

In HW equation of state, the relativistic electrons are 
assumed to be in /3-equilibrium with the nuclei (which 
also includes the neutron rich nuclei). Whereas, above the 
neutron drip density, the /3-equilibrium is among the nu- 
clei, free neutrons and the electrons. We have noticed that 
the presence of a strong quantizing back ground magnetic 
field modifies the /3-equilibrium configuration significantly 
although we have assumed that the electrons are only af- 
fected by the presence of strong magnetic field. As a con- 
sequence of which more heavier neutron rich nuclei are 
produced in presence of strong magnetic field in the inner 
crust region. In the conventional astrophysical scenario, in 
presence of free electron gas in stellar medium, the balance 
between Coulomb force and the nuclear force, which gives 
^^Fe as the most stable nucleus will get shifted towards 
the heavier nuclei. The nuclei, formed in such an environ- 
ment will contain more neutrons (by inverse /3-decay) com- 
pared to the usual picture. The Coulomb force plays a very 
little role. The nuclei becomes more and more neutron rich 



as the electron density goes up and when the matter den- 
sity increases to ~ 4 x lO^'^^g/cm'^, the ratio n/p reaches 
a critical level. Any further increase in density will lead 
to neutron drip in the medium. In such a scenario, highly 
neutron rich nuclei, electrons and free neutrons co-exist 
in chemical (/3) equilibrium. With the further increase in 
matter density, beyond neutron drip, nuclei, including the 
heavier ones melt and free neutrons appear in the medium 
and at some stage the kinetic pressure of the system will 
be dominated by the free neutron pressure instead of elec- 
tron pressure which may become negative if the magnetic 
field is strong enough. We shall see later that this posi- 
tive pressure contribution of free neutrons make the total 
pressure within the inner crust region a positive definite, 
even in presence of strong magnetic field. Whereas in the 
outer crust, the total pressure, coming from the electron 
gas only and is negative at the surface region of WS cells 
in presence of strong magnetic field. 

The BBP equation of state is applicable for nuclear 
matter in the density range from neutron drip density pdrip 
to normal nuclear density pnuc. {^ 2.8 x lO^^g cm^'^). It 
is well known that BBP equation of state is a consider- 
able improvement over the HW equation of state. In BBP 
equation of state, a mass formula, more or less like the 
HW equation of state is used but incorporated a lot of 
improvements from detailed many body calculation. The 
nuclear surface energy, for example, assumed in the pre- 
vious treatment to be that of a nucleus in vacuum. An 
introduction of free neutron gas out side the nuclei re- 
duces the nuclear surface energy. This is quite correct, 
because when inside and out side of a nuclei become iden- 
tical, the surface energy must vanish. In BBP equation of 
state, the nuclear Coulomb energy is included in a more 
accurate manner and is called nuclear lattice Coulomb en- 
ergy. It is well known that the BBP equation of state is 
applicable up to nuclear density, therefore, when all the 
bound nuclei dissolves into a continuous matter, mainly 
composed of free neutrons and a tiny fraction of protons 
and electrons, to incorporate this type of melting process 
of nuclei, in the BBP equation of state the factors giving 
the fractional volume occupied by the nuclei and the frac- 
tional volume occupied by the free neutron gas have been 
taken into account. It is also assumed that at this density 
(< Pnuc.) the nuclei are stable to /3-decay. Further, the 
neutrons in the free neutron gas are in chemical equilib- 
rium with the electrons and the nucleons within the nuclei 
and in addition to the /3-stability of the nuclei, the whole 
system must be in chemical equilibrium. In this model 
the kinetic pressure of free neutrons must necessarily be 
equal with the pressure of the nucleons bound within the 
nuclei. This is the condition for mechanical equilibrium. 
Since detailed mathematical formalism for both HW and 
BBP equation of states, which are assumed to be not af- 
fected by strong magnetic fields, are available in a large 
number of classic papers, including the original ones, and 
also included in many standard text books P^ITSIITB] . The 
effect of strong quantizing magnetic field on the electron 
part has already been discussed in the previous section. 
In this section, we therefore present only the numerical 
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results based on these two equation of states along with 
the results related to electron gas where ever needed, as 
developed in the first part of this section. 

In the numerical computation, we have found that for 
very low mass number, since there is no free neutrons 
available in the system, even much above the neutron drip 
density, only electrons contribute in total kinetic pressure 
and is negative if the magnetic field affects the electron 
part quantum mechanically. This result is almost identi- 
cal with the outer crust scenario. In fig. (7) we have plot- 
ted the variation of the critical value of mass number A 
and the corresponding atomic number Z with the strength 
of magnetic field using HW and BBP equation of states. 
These are the critical values at which the kinetic pressure 
of the inner crust region just becomes zero, i.e., the system 
just becomes mechanically stable. It is obvious from this 
figure that the minimum of mass number of the nuclei for 
which the inner crust matter just becomes mechanically 
stable, increases with the strength of magnetic field. In 
other words, this is equivalent to say that the presence 
of strong magnetic field makes the nuclei more and more 
massive or increases the minimum mass of the nuclei in 
the inner crust region. In some of our previous work, done 
long ago [12], we have seen that the same conclusion is 
applicable for quark matter. The presence of strong quan- 
tizing magnetic field generates quark mass dynamically in 
dense quark matter composed of massless quarks. This is 
the well known magnetic field induced chiral symmetry 
violation. The strong quantizing magnetic field acts like 
a catalyst to generate / increase mass. In the BBP equa- 
tion of state the minimum mass of the nuclei are more 
than that obtained from HW equation of state. Further, 
one should notice from this figure (figs. (7)) that for low 
and moderate strength of magnetic field, the effect is not 
so significant. But beyond lO^^G, when all the electrons 
occupy the zeroth Landau level, or in other words, when 
the quantum mechanical effect of strong magnetic field is 
most important, the minimum mass rises sharply to very 
large values. One can see from [T^] that the qualitative 
nature of the curve showing the dependence of dynami- 
cal quark mass on the magnetic field strength is exactly 
identical, although the physical scenarios are completely 
different. 

In fig. (8), we have plotted the variation of the ratio 
He/n and nn/n with the mass number using HW equa- 
tion of state. Here n,, is the electron density, Un is the free 
neutron density and n = n„ + n^A/Z is the total baryon 
density of inner crust matter. Dashed curve is for i? = 0, 
middle and the lower curves, indicated by el and e2 are 
for B — lO^Bc' and 10^i?c respectively. The curve indi- 
cated by n is for free neutron gas. In the HW equation of 
state, the neutron number density above the neutron drip 
point is independent of magnetic field strength. Further, 
the mass number around which neutrons are liberated 
from neutron rich nuclei is about 95, which is again inde- 
pendent of magnetic field strength. However, immediately 
after the emission of neutrons from heavy neutron rich nu- 
clei the overall kinetic pressure can not become non-zero. 
In fig. (9), we have plotted the same kind of variations as 



shown in fig. (8). Here we have used the BBP equation of 
state. Solid curves are for B = 0, while the dashed curves 
are for 10"^ Be . In this case the neutron drip out from 
heavy neutron rich nuclei around A = 100, which is lit- 
tle bit heavier and more neutron rich than the HW case. 
Further, the qualitative nature of x„ = n„/n is totally dif- 
ferent from HW case. Instead of increase initially and then 
saturates, as we observe in HW case, it decreases mono- 
tonically with A. Also, the free neutron density depends 
on the strength of magnetic field. However, the depen- 
dence is not so significant. These qualitative differences 
are found to be the consequence of chemical equilibrium 
among the free neutrons, electrons and nucleons within the 
heavy neutron rich nuclei and the overall /3-equilibrium 
condition. 

In fig. (10). the variation of total nuclear density is plot- 
ted against the mass number for both HW and BBP equa- 
tion of states. For the HW case, four different curves with 
magnetic field strengths: B = 0,B = 10b!^\ B = 10"^ B^"'' 
and B ~ lO^Bc are indicated in the diagram by 0, 1, 
10^ and 10^ respectively and shown by solid lines. In HW 
model, initially for low A values, only the bound nucleons 
within the nuclei (= An^/Z) contributes. Just beyond the 
critical value, A = 95, since neutrons drip out from the 
heavy neutron rich nuclei, the number density jumps sud- 
denly. Of course, below A — 60, we can not have stable 
inner crust matter. The same kind of variation for BBP 
equation of state are shown by the dashed curves and indi- 
cated by the symbols a, b and c. In this case, again because 
of chemical equilibrium conditions, the qualitative and the 
quantitative nature of the curves are completely different. 

The equation of state from both HW and BBP mass 
formulas are plotted in fig. (11). From these two figures it 
is quite obvious that in presence of strong quantizing mag- 
netic field, the inner crust matter becomes mechanically 
stable only at very high density, when enough number of 
free neutrons are available in the system. In this high den- 
sity situation positive neutron kinetic pressure dominates 
over negative electron pressure in presence of strong quan- 
tizing magnetic field. This is true for both HW and BBP 
equation of states. Further, it is also obvious from the fig- 
ures, that the matter becomes softer with the increase in 
magnetic field strength. The qualitative nature of both 
the equation of states are almost identical in presence of 
strong magnetic field. However, the HW equation of state 
is a bit softer than the BBP case. 



3 Conclusions 

In the first part of this article we have presented our in- 
vestigation on the properties of dense outer crust matter 
of the magnetars. We have replaced the dense metallic 
iron crystal by a regular array of spherically symmetric 
WS cells. It has been observed that the radius of each cell 
decreases with the magnetic field strength. In this article, 
however, we have not considered the magnetic field in- 
duced deformation of WS cells, which is quite likely if the 
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magnetic field is strong enough. In future we shall present 
this important issue with a completely different approach. 

It has also been noticed that the upper limit of Landau 
quantum number is a function of positional coordinate of 
the electron with which it is associated within the WS 
cells. We have observed that at the surface region, for 
all the values of magnetic field strength, this upper limit 
becomes identically zero. Which actually means that the 
electrons near the WS cell surface are strongly polarized in 
the opposite direction of external magnetic field. Whereas, 
for B > 10^^ G, they are polarized at every points within 
the cells. 

It has been observed that the electron density within 
the cells increases with the increase in magnetic field strength 
Further, for all the values of magnetic field strength, the 
electron density is maximum near the nuclear surface (r = 
r„) and minimum at the WS cell boundary (r = Tg). 

The electron kinetic pressure is found to be positive 
near the central portion of WS cell but is negative near 
the surface region. There is a position within the cell at 
which the kinetic pressure vanishes. The position of this 
point changes with the strength of magnetic field. We 
have also studied the variations of kinetic energy, electron- 
nucleus interaction energy, electron-electron direct poten- 
tial energy and electron-electron exchange interaction part 
within the cells. We have shown the variation of these 
quantities within the WS cells for a number of magnetic 
field strengths. 

In the second part of this work, we have investigated 
some of the properties of inner crust matter of magne- 
tars. We have used HW and BBP equation of states for 
the nuclear mass formula. It has been observed from both 
the models, that for a stable inner crust matter, the nuclei 
present must be heavier than iron and much more neutron 
rich. The heaviness is more in the case of BBP equation 
of state. We have also noticed that for low and moderate 
values of magnetic field strength, the variation of mass 
number and the corresponding atomic number with mag- 
netic field is not so significant. Whereas, for B > lO^^G, 
when electrons occupy only the zeroth Landau level, then 
much more heavier neutron rich nuclei are formed in the 
inner crust region. It is found that high magnetic field 
behaves like a catalyst to generates heavy neutron rich 
nuclei. 

We have observed that initially the electron density in- 
creases with the increase in mass number, but as soon as 
free neutrons appear in the system, the electron density 
decreases and saturates to some constant value which de- 
pends on the magnetic field strength. In the case of HW 
equation of state, free neutron density does not depend 
on the strength of magnetic field, whereas, for BBP case, 
because of chemical equilibrium condition, the free neu- 
tron density depends very weakly on the magnetic field 
strength. We have noticed that in the case of BBP equa- 
tion of state the overall qualitative difference is because 
of chemical equilibrium among the constituents. 

The total baryon density rises sharply like an avalanche 
for the value of A at which free neutrons appear in the 
system. However, for BBP equation of state, because of 



chemical equilibrium condition, the rise is not so sharply 
visible for a given magnetic field B. 

The qualitative nature of equation of states are almost 
identical. It is found that in presence of strong magnetic 
field, the inner crust matter becomes mechanically stable 
(with the positive value of kinetic pressure) only at very 
high density. 

We therefore believe that to study various properties 
of dense matter associated with magnetars, one must con- 
sider all these significant changes from the conventional 
scenario. Further, we expect that the strong magnetic 
field, if present well within the inner crust of magnetars, 
must affect the super-fiuidity of cold neutron matter. 
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Table 1. Variation of the parameter a and b for (j>{x) — ax + b 
and scaled surface radius Xs for WS cells, parameters a and 
/3 for Eke = ax'^ and the parameters p and q for E^^ = px'' 
with magnetic field strength 
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Fig. 1. The variation of surface radius fiXs in MeV ^ (solid curve) and the actual radius r^ in Aunit (dashed curve) with the 
strength of magnetic field (expressed in terms of Be ■) 
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Fig. 2. The variation of Vmux with x (dimensionless) . For the curves indicated by the symbols a and b, the relevant axes are at 
the bottom (x) and at the left side (i^max), whereas for the curves represented by the symbols a and /?, the axes are at the top 
and at the right side respectively. 
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Fig. 3. The variation of electron density in MeV"^ within the cells for B — h x Be with h — 1, 10, 100 and 1000 indicated 
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Fig. 4. The variation of —Een(x), the electron- nucleus interaction energy in MeV within the WS cells for B — h x Be with 
/i = 1, 100 and 1000. 
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Fig. 5. The variation of Eee{x), electron-electron direct interaction energy in MeV within the WS cells for B — h x Be with 
ft = 1, 100 and 1000. 
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Fig. 6. The magnetic field dependence of critical positions Xc (which is dimensionless) within the cell at which the total kinetic 
pressure vanishes, solid curve is for the numerically evaluated points and the dashed curve is the fitted functional form. 
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Fig. 7. Variation of critical values for A and Z with the strength of magnetic field B (expressed as B/Bc ) at which kinetic 
pressure becomes just zero. The solid curves are based on HW equation of state whereas the dashed curves are from BBP 
equation of state. In each case, the upper curve is for A and the lower one is for Z. 
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Fig. 8. Variation of the ratio nejn and nnjn with the mass number using HW equation of state. Here n^ is the electron density, 
n„ is the free neutron density and n = rin + neA/Z is the total baryon density of inner crust matter. Dashed curve is for 5 = 0, 
middle and the lower curves are for B 
for Un/n. 



,{e) 



j(e) 



10 Be and 10 Be indicated by el and e2 respectively. The curve indicated by n is 



18 Please give a shorter version with: \authorrunning and \titlerunning prior to \maketitle 

0.6 



0.4 - 



X 



<^ 



X 



0.2 



0.0 




I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I I 



60 



1 10 



A 



160 



210 



Fig. 9. Variation of the ratio Xe = rie/n and Xn = rin/n with the mass number A using BBP equation of state. Here rie is the 
electron density, n„ is the free neutron density and n = Un+neA/Z is the total baryon density of inner crust matter. Solid curve 
is for B — 0, the dashed curve is for 10'* x B^ . Electron part and the neutron part are indicated by Xe and x„ respectively. 
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Fig. 10. Variation of mass number of the nuclei with the total baryon density. HW equation of state: four different magnetic 
field strengths:^ = 0, B = loM"', B = lO'^Bi"^ and B = 10^5^^', indicated by 0, 1, 10^ and 10^ respectively. BBP equation of 



state: Dashed curves, the magnetic field strengths B = 0, B = lOBc (both are indicated by a), _B = 10 Be (indicated by b) 
and B = 10'^ B^ (indicated by c). Here, the lower x-axis is for HW equation of state and the upper one is for BBP equation of 
state. 
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Fig. 11. Equation of states from HW and BBP mass fornmlae. Sohd curves arc for HW equation of state, whereas, the dashed 
curves are for BBP case. For each case, the upper curve is for B = 0, middle and lower curves are for B — Be and B — 10'^ Be 
respectively. Please note that the scales for a;-axes at the top and bottom are different because of two different equation of 
states. 



